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We study linearized equations of a ghost-free gravity in four- and higher-dimensional spacetimes. 
We consider versions of such a theory where the nonlocal modification of the □ operator has the form 
□ exp[(—where V = 1 or V = 2n. We first obtain the Newtonian gravitational potential 
for a point mass for such models and demonstrate that it is hnite and regular in any number of 
spatial dimensions d > 3. The second result of the paper is calculation of the gravitational field 
of an ultrarelativistic particle in such theories. And finally, we study a head-on collision of two 
ultrarelativistic particles. We formulated conditions of the apparent horizon formation and showed 
that there exists a mass gap for mini-black-hole production in the ghost-free theory of gravity. In 
the case when the center-of-mass energy is sufficient for the formation of the apparent horizon, the 
latter has two branches, the outer and the inner ones. When the energy increases the outer horizon 
tends to the Schwarzschild-Tangherlini limit, while the inner horizon becomes closer to r = 0. 

PACS numbers: 04.70.-s, 04.50.+h, 04.50.Kd 


I. INTRODUCTION 

Singularities are inherent properties of general relativ¬ 
ity. It is generally believed that the Einstein-Hilbert ac¬ 
tion should be modified in spacetime domains where the 
curvature becomes large. Such a modification is required, 
for example, when one includes in the theory quantum 
corrections, connected with particle creation and vacuum 
polarization effects. At a more fundamental level, the 
modification of the gravity equation might be required if 
the gravity is described as an emergent phenomenon. In 
such a case the Einstein equations are nothing but the 
low energy limit of the corresponding more fundamental 
background theory. The string theory is a well-known 
example. It is convenient to introduce two (generally 
different) energy scale parameters /i and ft. The corre¬ 
sponding length scales are A = and A = . We 

assume that when the spacetime curvature fZ is much less 
than A“^, the corrections to the Einstein equations are 
small. These corrections become comparable with other 
terms of the Einstein equations at 7^ ^ A“^, and for 
higher values of the curvature they play an important 
role. We assume that one can use the classical metric 
for the description of the gravitational field. Eor ex¬ 
ample, one can understand it as a quantum average of 
some metric operator, = {g^y). This means that the 
quantum gravity effects, and in particular fluctuation of 
the metric, are small. In other words, one can use the 
effective action approach to study spacetime properties 
in this domain. The second parameter. A, defines the 
scale when effective action description breaks down and 
the quantum nature of the gravitational field becomes 
important. 

In studies of the singularity problem in modified grav¬ 
ity it is usually assumed that fi ^ g. In the present 
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paper we also use this assumption and discuss some as¬ 
pects of the singularity problem in the framework of the 
classical modified gravity equations. 

There exist a wide class of the modified theories of 
gravity proposed to solve fundamental problems of black 
holes and cosmology. We consider a special class of such 
theories, namely theories with higher derivatives. Impor¬ 
tant features of such theories can be clarified already in 
a simple approximation when the gravitational field is 
weak and can be described as the perturbation on the 
fiat spacetime background. Such an analysis was per¬ 
formed by Stelle [T]. In particular, he demonstrated that 
the Newtonian gravitational potential of a point mass lo¬ 
cated at r = 0 can be made finite at this point, if the 
higher derivative terms are included in the gravity equa¬ 
tions. Detailed analyses of this problem can be found in 
recent papers El [31. 

However, the higher derivative gravity, as well as any 
theory with higher derivatives, has a fundamental prob¬ 
lem. In a general case the propagator of such a theory 
contains two or more poles, and, as a result, it almost 
always contains ghostlike excitations (see, e.g., [T] and 
mM)- Presence of the excitations with negative energy 
results in an instability of the theory and the possibility 
of an empty space decay. This is a special case of a very 
general phenomenon known as Ostrogradsky instability 
[6] (see discussion in [5]). 

In the higher-derivative theory a standard box opera¬ 
tor □ , which enters the field equations is changed to the 
operator P(n)n, where P{z) is a polynomial. The poles 
of P~^{z) correspond to additional degrees of freedom. 
However, there exists an interesting option of theories 
where P~^{z) is an entire function of 2 ; and hence it does 
not have poles in the complex plane. Such a modification 
of the gravitational equations is called ghost-free (GE) 
gravity (see, e.g., IMl] and references therein). GE grav¬ 
ity contains an infinite number of derivatives and, hence, 
it is nonlocal. Theories of this type were considered a 
long time ago (see, e.g., [T6H2D]). They appear natu- 
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rally also in the context of noncommutative geometry 
deformation of the Einstein gravity mi ns (see a review 
[23] and references therein). The initial value problem 
in nonlocal theories was studied in •The applica¬ 

tion of the ghost-free theory of gravity to the problem of 
singularities in cosmology and black holes can be found 
in [25143T] . Static and dynamical solutions of the lin¬ 
earized equations of the ghost-free gravity in four and 
higher dimensions were studied in [32l [33] . Recently the 
consequences of the ghost-free modifications of higher¬ 
dimensional gravity on the entropy of black holes and on 
cosmological models have been studied [34|. 

In this paper we continue study of the linearized equa¬ 
tions of the GF gravity. In Secs. [npV| we study solutions 
for a static gravitational field in the Newtonian approx¬ 
imation in different models of the GF gravity. Namely, 
we consider a class of the GFn theories of gravity with 
P(n) = exp[(—A static solution of the lin¬ 
earized equations for A" = 1 in four-dimensional space- 
time was found in [8j[26] (see also [33|). In this paper we 
generalize this result to the higher-dimensional case and 
obtain new solutions for GF 2 n theories in the spacetime 
with an arbitrary number d of spatial dimensions. In Sec. 
[V|we used these results to obtain a solution of the GF 
gravity describing a gravitational field of an ultrarela- 
tivistic particle. We succeeded to find a generalization of 
the famous Aichelburg-Sexl solution [35] to the GF grav¬ 
ity in an arbitrary number of dimensions. In Sec. jVIj we 
used the obtained solutions to study the apparent horizon 
formation in head-on collision of two ultrarelativistic par¬ 
ticles. This problem for the general theory of relativity in 
four dimensions was first solved by Penrose [36]. Later, 
this result was generalized for a collision with a nonzero 
impact parameter in four and higher dimensions [STHlO] . 
In the present paper we show that in the GF gravity a 
similar process has two important new features: (i) the 
apparent horizon is not formed if the center-of-mass en¬ 
ergy of the particles, P, is smaller than some critical 
value Ecrit^ which depends on the scale parameter /i, the 
type of the theory, and the number of spacetime dimen¬ 
sions; (ii) if the energy is larger than Ecrit the apparent 
horizon besides the usual outer part always has another 
inner branch. We discuss the obtained results in the last 
section. 

In the present paper we use units in which h = c = 1 
and sign conventions adopted in the book mi¬ 


ll. NEWTONIAN LIMIT OF 
HIGHER-DIMENSIONAL HIGHER-DERIVATIVE 
EQUATIONS 


Let us consider a static gravitational field perturbation 
on a flat background and write the corresponding metric 
in the form 


ds^ = —(1 -h 2(p) dt^ -h (1 — 2-0 -h 2ip) di^^ 
di^ = 6ik dx^dx^ , X* = {x^ ,..., . 


Here and later we denote by d = D — 1 a number of 
spatial dimensions. We also have 


h/QO — ^ij — ‘2‘ijP 

h = 2[(d -h l)(f — d^p] . 


By su bstit uting these expressions into the gravity equa¬ 
tions (A5) one gets 

a(A) A'ip = Kd{Too + 

[a(A) - (ic(A)] Aip+ {d- l)c(A) A t/) = HdToo- 


Here Kd = and D = d + 1 is the total number of 

spacetime dimensions. In the Newtonian approximation 
S'^^Tij = 0 and the first of these equations takes the form 

a(A) A^p = KdToo . (2.2) 


For the gravity theory with c = a the equations simplify 
and one obtains 


d-l 


(2.3) 


and the metric (2.1) takes the form 


ds^ = —(1 -h 2ip) dt^ + (1 — ^—^^) • (2.4) 


For a point mass m the energy density has the form Too = 
m5^(x). Then for the Einstein gravity, where a = c = 1 
one has 


ip = 


Kcjmr (f) 


1 


2(d — l)7r^/^ ^ 

In four dimensions D = 4 (d = 3) 

hzs m 


ip = - 


Stt r 


(2.5) 


( 2 . 6 ) 


III. STATIC SOLUTIONS OF LINEARIZED 
EQUATIONS IN GHOST-FREE GRAVITY 

A. Ghost-free gravity 


The Newtonian potential (2.5) is evidently singular at 


r = 0. One can regularize it and make it finite at r = 0 by 
modifying the gravity equations in the ultraviolet (UV) 
domain. For example, one may assume that a(n) and 
c(n) are polynomials of the □ operator. If these func¬ 
tions obey the condition a(0) = c(0) = 1, the theory 
correctly reproduces the standard results of general rela¬ 
tivity in the infrared regime, that is in the domain where 
r ^ oo. In a general case such a theory possesses ghosts. 
These ghosts are new degrees of freedom which are con¬ 
nected with extra poles of the operators a~^ and c 
which give contributions to the propagator with a wrong 
(negative) sign. However, there exists an option to use 
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such functions a~^{z) and c~^(z) that are entire func¬ 
tions of the complex 2 ;-variable which do not have poles. 
It happens, for example, when a{z) and c{z) are of the 
form exp(P(^)), where P{z) is a polynomial. A modified 
gravity which contains such regular formfactors is called 
ghost-free (GF) gravity. In the present paper we focus on 
the special class of the theories of GF gravity. Namely, 
we assume that 

a{n) = c(n) = expii-n/fiY) • (3.1) 

We denote such a theory GF^. We restrict ourselves by 
considering the cases N = 1 and N = 2n, which are of 
the most interest for applications. 

The exponent of the operator can be written in the 
form of a convergent series of the powers of this opera¬ 
tor. However, it is not a good idea to “approximate” the 
exponent by the polynomial which is obtained by keeping 
a finite number of terms in this series. The inverse opera¬ 
tor will have extra poles and the ghost will be present for 
such truncation. That is why our first goal is to present 
these nonlocal objects in the form of an integral trans¬ 
form which contains a well-defined kernel. 


B. Potential fjd and Green functions in GF theories 


Gonsider the equation for the potential fjd created by 
a point massive particle placed at a point x' 

Ftpd = Heims’^{x - x'), (3.2) 


where the operator F is defined on the d-dimensional 
Euclidean space. It is assumed to be a function of the 
Laplace operator 

F = F{-A), FiO = -ea(-e- (3.3) 

The Euclidean Green function Dd{x,x') of this operator 
is the solution of the problem 

FDd{x,x') = -6^{x-F) (3.4) 


with vanishing boundary conditions at infinity. Eormally 
it can be treated as a matrix element 

Dd{x,x') = {x\D\x') (3.5) 


of the operator 

D = 


D = D{-A), 


F{0 Cai-0' 


(3.6) 


The momentum space calculations of D djx^ x') are 
seated in Appendix]^ The result reads ( |B6[ ) 

X Jd_dVv\x - x'\), 


pre- 


(3.7) 


In Sec. El and Sec. EH we will use this Green function 
to study a gravitational field created by ultrarelativistic 
particles. Eor this purpose it is useful to have another 
representation of the Green function, where the Bessel 
function is replaced by its integral representation 


1 rc+ioo 

■’-(‘I = (2) ssL,. 


c—zoo 

c > 0. 


expit-j 


Then after the change of the integration variable 

t = irjr, > 0, (3.8) 

the Green function can be written in the form 


1 ~ 

Dd{x,x') = —J driD{ri) 

/ oo—ic 

. 00 -ic (47rir)<^/2 




(3.9) 


Note that the last integral contains the expression which 
is known as the heat kernel of the Laplace operator in a 


d-dimensional fiat Euclidean space 


f /IN 1 

Kd{x,x r - . e 4 . . 

(3.10) 

The heat kernel obeys the equation 


idrKd{x,x'\T) + AKd{x,x'\T) = 0 

(3.11) 

and the condition 


lim Kd{x, = ^^(^ “ ^ 0 * 

r^O 

(3.12) 

It describes the amplitude 


Kd{x,x'\T) = {x\F'^^\x'). 

(3.13) 


In flat space, because of the symmetries of the system 
in question, both the Green function Dd and the potential 
fjd are the functions of a distance r between the points 
only 


Dd = Dd{r), tpd = ‘>Pd{r), r = ^/(x - x'Y. (3.14) 

The potential at the point x created by the massive par¬ 
ticle located at the point x' is 


f)d = -ndmDd{r). (3.15) 

IV. GRAVITATIONAL POTENTIAL IN 
LINEARIZED GF GRAVITY THEORIES 

A. General properties of GF theories 

All GF theories of gravity are assumed to reproduce 
Einstein gravity in the low energy regime, i.e., at large 
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scales. In particular it means that the functions a(^) and 
c(0 approach smoothly to 1 at small 

a(0 = 1 + 0(0, c(0 = 1 + 0(0. (4.1) 

Then we have the functions F{^) = —^+0(0 and D{^) = 
1/.^ + 0(1). This property and (3.7),(3.15) guarantee 


that in the limit of large distances one gets a universal 
asymptotic for the potential for all these GF theories: 


ipdir) 


r(f-i) 


(4.2) 


Obviously, as it sho uld b e, it exactly reproduces the grav¬ 
itational potential (2.5) in the higher-dimensional Ein¬ 
stein gravity theory. 

The asymptotic of the potential at small distances is 
theory dependent. Our particular interest is in GF^ the¬ 
ories, where 


a(-e)=exp((e/M^)^) 


(4.3) 


and A/" = 1 or an even integer number. The parameter 
fi characterizes the scale where the nonlocality becomes 
important. One can show that for all GF^ gravities the 
potential ipd is finite at small r. For these theories the 
asymptotic at rji 0 can be compute d explicitly. Let 
us substitute ( |4.3| ) to ( |3.6[ ),(3.7),(3.15) and change the 
integration variable rj = z^/r+Then we have 


'ipdir) = - 


Kdm 


(27r)^/^7 


.d-2 


f 


dzz^-^e Jd_i(z). 


(4.4) 

One can see that in the limit when r/i ^ 0 only small 
argu ment s of the Bessel function contribute to the inte¬ 
gral (4.4). Therefore, one can substitute there an expan¬ 
sion 


Jd_i(^) = 


(I)- 


nd/2) 


i-S 


2d Sd{d + 2) 


+ 0{z^) 


Then taking the integrals in (|4.4|) one obtains 


(4.5) 


,d-2 


ipdir) - Kam- 


Y (^^) - (—) 

^ \ 2N ) 2d ^ \2N) 


(47r)+27vr (f) 


-O(rV). 

(4.6) 

One can see that the leading term is finite and propor¬ 
tional to . Moreover the next term in the ex¬ 

pansion is proportional to that guarantees regularity 
of the metric at r = 0. 

There are other interesting universal properties of the 
potentials in generic GF graviti es. F or example, because 
the distance r in the integral (3.7) does not enter the 


function D and due to the properties of the derivatives 
of Bessel functions it is clear that there is a universal 
relation 




(4.7) 


For the potentials, considered as functions of the radial 
distance r, this property leads to the relation 

-J—i’d+2{r) = —(4.8) 
Kd +2 Kd 27rr dr 

provided the mass parameter m is the same in d and 
{d + 2) dimensions. 


B. Potential in GFi theory 

The static potential ipd in the GFi theory satisfies the 
equation 

exp(— A /A'lpd = K,d'm6^{x — x'), (4.9) 

so that 


F{0 = = 






(4.10) 


Substitution of this expression into (3.7) and change of 
the integration variable r] — z^ leads to 


Fd{T) = 


(27r)^/^r 


2rr>d—2 


[ 

/^j^d/2rpd 2 


dzzi-^e~^Ji_dz) 


(4.11) 

where ■j{n,x) is the lower incomplete gamma function 
m- At large distance r ^ ji ^ this expression re¬ 
produces the static Green function of the d-dimensional 
Laplace operator 


Gd(x x') = ^ _F 

'^d(X,X ) ^^d/2j.d-2' 


(4.12) 


For small distances r ji ^ the Green function Dd{r) is 
a regular function of r and is of the form 

Dd{r) = ,, ^C.\d^2 (1 - +■■■ ■ (4.13) 


{d — 2)(47r)^/2 V ^ 

The potential ipd is given by 

^d = -KdmDd{x,x') = -Kdm ■ (4.14) 

In four-dimensional spacetime (d = 3) we reproduce the 
results of [8l|23l|26l®ll3] 

erf (r/u/2) 


V’3 = -Ksm ■ 


Attt 


(4.15) 


In the case of five-dimensional spacetime (d = 4) we ob¬ 
tain even simpler expression 

1 — exp (—r^/i^/ 4 ) 

n = -. (4.16) 

The potentials jjjd in an arbitrary number of dimen¬ 
sions qualitatively look alike. They are negative and fi¬ 
nite at r = 0. At larger distances they become more 
shallow and at r ^ quickly approach the Einstein 
asymptotic (4.2). 
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C. Potential in GF2 theory 


When N = 2 the operator F corresponds to 

a(A) = exp(A^//i^) (4.17) 

and, hence, 

F{i) = (4.18) 

Then the potential takes the form 


i^dir) = - 


d{d- 2 ) 2 ^-^n^ 


X 



/d 1 1 d d 
V4~ 2’2’4’4^ 



^id- 2 )y (d 3 d d,l.\ 
r(|+i)^ T4’2’4^ ’4^2’^ )\ ’ 

(4.19) 


where 


y = 


r'^lb 


16 


(4.20) 


and pFq is the generalized hypergeometric function (see, 
e.g.,|42]). 

Qualitatively the potentials for different parameters N 
and in different dimensions d look similar. Figs. [T][^ show 
examples of the gravitational potential for d = 3 and 
d = 4 in two cases, N = 1 and N = 2. 


theory. For all these theories the asymptotic at large 
distances is governed by the (4.2) and the asymptotic at 
small distances is given by (|r^ Let us present here only 
one more explicit example of the potential in GF 4 gravity 


d(d + 2)(d + 4)2r 

m 

iBi 

- 8 y{d + 2 ){d + 4)2r 

1 B 2 

+32/(d + 4)^r 

00 —[— 

to 

1 Bs 

2048 3 ^ / 

'd+ 12^ 


- 3 


1 B 4 


(4.21) 


where 


= F fd^.l 1 3 d d|3 dll d I 1 . \ 

l-^7\^ 8 ’4’2’4’8’8“'“4’8“'“2’8“'“4’ 256 J ’ 

Bo = ^ ^ ^ + ^ ^ + ^ d + 

J-'Z 1-t 7\^8’2’4’4’8^-^’8^ 4’8^ 2’8^ 4’ 256 J ’ 

Ro = F f^+2.3 5 3 d 1-1 d I 3 dll d , 5.R\ 

1-^ 7 8 ’ 4 ’ 4 ’ 2 ’ 8 8 4 ’ 8 2 ’ 8 4 ’ 256 y ’ 

Ba = F 3 7 d i-i d I 3 d I 3 d I 5.R\ 

-^4 8 ’ 4 ’ 2 ’ 4 ’ 8 “T-L, 84 5 82 ’ 8 4 1 256 y ’ 

and the coefficient 


“ 22<i-57r^ d{d - 2){d + 2){d + 4) 

1 

"GW(i)F(wm‘ 

Expressions for the potentials become more complicated 
for higher N and we do not present them here. 


D. Potential in GFn theories 

Similar results in terms of the generalized hypergeo¬ 
metric functions can be derived for an arbitrary GFn 


V. PENROSE LIMIT 

Let us demonstrate now, that obtained static solutions 
of the GF gravity can be used to find the gravitational 
field of an ultrarelativistic object. In the standard 4D 
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Einstein gravity such a limiting metric is known as an 
Aichelburg-Sexl metric [35] . This metric was generalized 
to the case of higher dimensions and for the spinning ob¬ 
jects (called gyratons) in papers |44]-[46]. In this section 
we obtain a metric created by an ultrarelativistic object 
moving in T)-dimensional spacetime (nonspinning gyra- 
ton metric) in GF theories of gravity. As we shall see 
a key role in th is de rivation is played by the heat kernel 
representation (3.9) of the Green function Dd{x^x'). 


Consider the metric in the following form 

ds^ = -(1 + 2tfd) dt^ + (1 - 2'tpci + 2ipd){dy^ 

x = {y,U), a = 


del), 

(5.1) 


Let us boost this metric in the ^-direction 

i=l{y-l3t), y = j{t-l3y), 7 = (1 -(5.2) 

and introduce null coordinates 

u = t-y, v = t-\-y. (5.3) 

In the relativistic limit, when the boost velocity is close to 
the speed of light, i.e., P ^ the boost factor 7 ^ oo. 
In this li mit dt ^ ^du and dy ^ —^du. Then the line 
element (5.1) becomes 

ds^ = —dudv + d(^\ + ^ (5.4) 

where 


= -2 lim (y^V'd)- 

7^00 


(5.5) 


For a point particle of mass m the Penrose limit corre¬ 
sponds to ultrarelativistic limit 7 ^ oc with the condi¬ 
tion that an energy E = 7 m of the pa rticle i s ke pt fixed. 

The gravitational potential 'ipd (see ( |3.15 ),( [3^ ) can be 
presented in the form 


V’d = -KdmDdir) 

dyD{y) j 


dr 


iF (5.6) 


27r 


(47riT)'^/2 


One can see that the boost affects only the last exponent 
in this integral representation. 

Taking into account that after the boost 


r2^7'(W)'v, y = (c±-c;)t 

and using the delta-function representation 


lim ^ ^ = S(u), 

7^00 7/47rir 

= Fd(p)S(u - u'). 


we obtain 


Here 


Fd{p) = '^ dyb{y) 


dr 


(dTrir)^^ 


(5.7) 


(5.8) 


(5.9) 




(5.10) 


Comparison of this integral expression with (5.6) leads to 
the observation that the function Fd{p) is proportional 
to the gravitational potential defined in space of one di¬ 
mension less, i.e., in the space orthogonal to the particle 
motion: 


d^d{p) = ‘^f^dE Dd-i{p) = -2 


f^dE 

Kd-im 


i^d-i{p)- 


(5.11) 


This property is valid for arbitrary GF^ theories of grav¬ 
ity. 

Using the property (4.8), which is also valid for a 
generic GF^ gravity, we derive a relation 


d,Ed{p) = ArEpijdriip). E = (5.12) 

Kd+im 

This relation will be useful for the study of gravitational 
effects in collisions of ultrarelativistic particles (gyratons 
|45H47|) in the next sections. 


VI. APPARENT HORIZON FORMATION FOR 
HEAD-ON COLLISION OF THE 
ULTRARELATIVISTIC PARTICLES 

Our next goal is to use the obtained results to study 
head-on collision of the ultrarelativistic particles in the 
GF theories of gravity. We use an approach developed by 
Penrose [36] and D’Eath and Payne [48ll5Q] and approx¬ 
imate the colliding particles by gyratons. A schematic 
picture of such a process is shown in Fig. [^ It shows 
two-particle motion in the center-of-mass frame. Each of 
the particles moves with the velocity of light. Particle 
1 moves from the left to the right along the ^-direction, 
while particle 2 moves in the opposite direction. The null 
lines, representing their trajectories, belong to = 0 and 

V = 0 null planes, correspondingly. The gravitational 
field of these particles is localized on the plane = 0 (for 
particle 1) and 1 ; = 0 (for particle 2). The intersection 
of two null planes is the {d — 1 )-dimensional transverse 
plane. In the regions /, //, and ///, outside the u = 0 
and 'T’ = 0 null planes the metric is flat and null rays in 
these domains are nothing but null straight lines. How¬ 
ever, when such a ray passes either through rt = 0 or 

V = 0 planes, it is scattered by the gravitational field of 
the corresponding particle. 

Our purpose is to study formation of the apparent hori¬ 
zon in such a process. Let us remember that a trapped 
surface is a compact spacelike {d — l)-dimensional sur¬ 
face which has the property that both of the null con¬ 
gruences orthogonal to it, are not expanding. We focus 
on the outgoing congruence. One calls a trapped surface 
a marginally trapped surface if the outer normals to it 
have zero convergence m- In a spherically symmetric 
spacetime one may consider spherical slices and define 
an apparent horizon as d-dimensional surface which on 
each of the slices coincides with the marginally trapped 
surface. 

















7 



FIG. 3. Head-on collision of two ultrarelativistic particles. 


The problem of ultrarelativistic particle collision in 
general relativity was discussed recently in connection 
with possible mini-black-hole creation in colliders [S3- 
EgiiT]. Eardley and Giddings [37] demonstrated that a 
problem of existence of the apparent horizon can be re¬ 
duced to a special boundary-value problem for an elliptic 
(Poisson) equation in a flat spacetime. Generalizations of 
these results to the collision of shock waves on AdS back¬ 
ground were also considered in [52] [53] . The problem is 
greatly simplified for the case of the head-on collision and 
can be solved analytically in any number of spacetime di¬ 
mensions. In the present paper we follow their approach. 
Let us write the metric ( |5.4| ) in the form 

ds^ = —dudv + d^\ + ^ddu^ , ( 6 . 1 ) 

$ = P = (6.2) 

It is possible to show that geodesics and their tangent 
vectors are not continuous in these coordinates (see e.g. 
m)- One can change the coordinates so that both 
geodesics and their tangent vectors will be continuous in 
the new coordinates. The new coordinates in the domain 
II are defined as follows 

— u 

U = U, Ci = Ci + , 

1 (6.3) 

V = V + $'d{u) + -ui9(m) (V^?)^ . 

A similar transformation (with a change u v) should 
be made in th e do main III. 

The metric ( |6.1| ) in the new coordinates takes the form 

ds^ = -dudv + - 5i^]dCidCk, 

(6.4) 

We consider a special marginally trapped surface S 
which consists of two parts Su Sy. In coordinates (u^vXi) 


a position of Su and Sy on two incoming null planes is 
described by equations 

{u = —T(p), = 0} and {i 4 = —T(p), u = 0} , (6.5) 

respectively. These two {d — l)-dimensional surfaces in¬ 
tersect at {d — 2 )-dimensional boundary C, located at 
u = V = 0. The function T is positive inside the bound¬ 
ary C and vanishes at C. The internal (induced) geometry 
of Su and Sy are the geometry of a half of a (d — 1 )- 
dimensional round sphere, their intersection C being a 
round {d — 2)-dimensional sphere. For the head-on col¬ 
lision the function T(p), which enters both equations in 
( |6.5| ), is the same. In m it was shown that the outer 
null normals have zero convergence in Su and Sy if 


V2(«'-J-d)=0. (6.6) 

A condition that both normals (in Su and 5^) coincide 

at their boundary C implies 

(V«')2=4. (6.7) 

Denote y = T — and by pc the radius p at the 
boundary. Then 

vy = 0 , xc = -:Fd{pc) . ( 6 . 8 ) 

Hence one can put X = ~d^d{pc) inside C so that 


^ = J^d{p) - J^d{pc) • 
The condition ( |6.7| ) takes the form 
(VJd)%=4. 


Using (5.12) one gets 

2'KKdE pDd+1 (p) = I- 


(6.9) 


( 6 . 10 ) 


( 6 . 11 ) 


In terms of a dimensionless coordinate x = /ip, dimen¬ 
sionless energy E = 27t p,^~‘^ K,d and a dimensionless 
profile function 


Pd{x) = xDd^i{xlp)lp^ 


d-l 


this condition reads 


( 6 . 12 ) 


Pd{x) = 

h/ 


(6.13) 


All functions Pd{p) look similar (see Figs.an d§. They 
vanish at x = 0 and then grow, reach maximum, and then 
decrease to a universal asymptotic, that does not depend 
on the parameter though depen ds on d. The plots 
Figs. and show solutions of Eq.(6.13). The appar¬ 
ent horizon exists for the energy obeying the condition 
E > E^^mcai- In this energy domain it has at least two 
branches, inner and outer. At E = they meet and 

the apparent horizon disappears [54]. This behavior re¬ 
sembles qualitatively that of the colliding relativistic ex¬ 
tended sources [55]. This resemblance is not accidental. 




PsM P4M 




X 


FIG. 4. The plot shows function P^{x) for N=l,2. 


FIG. 6. The plot shows function Pa{x) for N=l,2. 


X 



FIG. 5. The plot shows the radius x = /xp of outer and inner 
apparent horizons as a function of the energy E for d = 3 and 
N=l,2. 





FIG. 7. The plot shows the radius x — jip oi outer and inner 
apparent horizons as a function of the energy E for d = 4 and 
N=l,2. 


One can rearrange Laplace opeators in (3.2) and move 
a(A)“^ to the right-hand side of the equation. Then it 
can be identically rewritten as 

At/>d=j, j = Kdma{A)~'^S'^{x-x'), (6.14) 

When acting on the localized source, the operator 
a{A)~^ delocalizes it and makes j to become effectively 
an ex tended current for the traditional Laplace equation 
(6.14). In this sense the analogy of effects in the ghost- 


free gravities and for the colliding extended sources m 
becomes evident. 


VII. SUMMARY AND DISCUSSION 

In this paper we discussed an application of the lin¬ 
earized equations of the ghost-free theory of gravity to 
three connected problems. First, we calculated the grav¬ 
itational potential of a point mass in the Newtonian limit 
and showed that GF modification of gravity works as a 
regularizer. Namely, this potential is regular at the ori¬ 
gin. This property is valid for GFi and GF 2 n theories 
in any number of spatial dimensions d > 3. This is a 
generalization of the earlier obtained result for GFi for 


d = 3 [H [26] and for d > 3 [33| . The second main re¬ 
sult of the paper is calculation of the gravitational field 
of an ultrarelativistic particle in the GF^ theories. The 
obtained metrics are generalizations of the famous four¬ 
dimensional Aichelburg-Sexl metric m of general rela¬ 
tivity. Again, the obtained metrics are solutions of the 
equations of the GF^ gravity equations {N = 1 and 
N = 2n) in a spacetime with an arbitrary number of 
dimensions d + 1. And finally, we used these results to 
study an apparent horizon formation in the head-on col¬ 
lision of two ultrarelativistic particles. Our main conclu¬ 
sion is that in such a process there exists a mass gap for 
the mini-black-hole formation. If fi is the characteristic 
mass scale of the corresponding ghost-free theory, then in 
order for a mini-black hole to be formed in the collision, 
the center-of-mass energy E should be of the order of or 
larger than . Another important feature of 

the process is that when the apparent horizon is formed, 
it has two branches: outer and inner marginally trapped 
surfaces. Both of them have the geometry of the sphere. 
When the center-of-mass energy increases, the inner part 
becomes closer to the point until it reaches the scale A, 
where the model we used breaks down. 

This result is again valid for any GF^ theory {N = 1 
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and N = 2n) in any number of dimensions. It can be con¬ 
sidered as some indication that for such theories the in¬ 
ner singularity of a black hole might be absent and there 
exists a closed apparent horizon. Such a model was pro¬ 
posed in [56] and discussed later in many publications. It 
should be emphasized that most of the results, related to 
the study of the models with closed apparent horizons, 
beyond a linear approximation, were obtained without 
using concrete dynamical equations. In this sense they 
are phenomenological. It is a real challenge to obtain so¬ 
lutions for a dynamical collapse in the modifications of 
the Einstein theory which are UV complete. In particu¬ 
lar, in order to arrive at a definite conclusion concerning 
the structure of a black hole interior in the GF gravity 
one needs to perform analysis in the complete version of 
such a theory, which includes nonlinear effects. 
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Appendix A: Linearized equations of the 
higher-derivative modification of the gravitational 
equations in higher dimensions 


In order to obtain linearized equations of a theory of 
gravity with higher derivatives in higher dimensions one 
can follow a similar derivation in four dimensions pre¬ 
sented in the papers laiia. In this appendix we collected 
the corresponding formulas for further reference. 

The main steps of this derivation are the following. 
One considers first a covariant action which besides the 
Einstein term contains also a part Sq which is quadratic 
in curvature. The latter may contain an arbitrary num¬ 
ber of covariant derivatives acting on each of the curva¬ 
ture tensors. One can always move the derivatives acting 
on the first Riemann tensor to the position, when it acts 
on the other one. This can be achieved by using integra¬ 
tion by parts. The number of derivatives may even be 
infinite, so that a theory is nonlocal. Since each of the 
Riemann curvature tensors has fourindices, the maximal 
total number of derivatives with “free” indices is eight. 
All other derivatives can be combined in functions of the 
covariant box operator. In order to achieve this it might 
be required to commute the derivatives. But this opera¬ 
tion produces terms which are of the third order in the 
curvature so that they should be neglected in the adopted 
approximation. Using symmetry properties of the curva¬ 
ture tensor, Bianchi identities and commutativity of the 
covariant derivatives in the adopted approximation one 


finally obtains the following expression for Sq [Siiig 
5 = dx^[R + RFi{n)R + R^^F2{n)R^'' 

+R^uXaF3{D)R^^-'^-] . 

Here Kd = Sir and is the gravitational coupling 
constant in D-dimensional spacetime. In four dimensions 
the value of this constant is fixed by the requirement that 
the Poisson equation for the gravitational potential in the 
Newtonian limit has a standard form. There is an ambi¬ 
guity in the normalization of G^^^ in higher dimensions. 
We fix it by requiring the Einstein-Hilbert action to have 
the same form in all dimensions. 

This general form of the quadratic in curvature action 
can be further simplified using the following observation 
ElEZ]: the ” Gauss-Bonnet structures” of the form {k > 
1 ) 


+ i?tfi? = 0(i?3) + div. 

(Al) 

in arbitrary dimensions are all of the third and higher or¬ 
der in curvature plus total divergence terms. As a result, 
the general higher derivative action can be written in the 
form which contains only two arbitrary functions of the 
box operator [12]. 

To obtain the linearized equation we write the action 
in the form 

S=^{So + Si + S2 + Si + S 3 ), 

2nd 

^0 = y dx^/^ R, 

Si = J dxy^RFi{a)R, 

82 = J dx^R^,F 2 {a)Rf^‘', 

83 = J dxV^R^,x„F3{a)R'^''^F 

We use the following expressions for the variations of the 
objects that enter the above action and keep only the 
terms that are quadratic in perturbations 


/' 


So = - I dx{ h° 


1 , 


+ -hUh , 


Si = J dx Fi{U)d^dM& 

h + hn^Fi{n) h ), 
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S2 = \J dx{2h^^F2{a)d^dMph"^ 

-2h^^'' □F2(n)a^a„ - 2h^^^ UF2{U)d^d^ h 

+h^'' □2^2(0) V + h o^F2in) h), 


for static gravitational potentials, which contain contrac¬ 
tions of the form r]^yh^^^ would be explicitly dependent 
on D. 

Appendix B: Gravitational potential in momentum 
space 


Ss = Jdx F3{n)d^dM0 

□2F3(n) nF3{n)d^da h'^u) ■ 

Let us write the total linearized action S in the form 


— f dx(- 

2 KdJ V 2 


S=— I dxi aahf,^ + hf^’^bd^da h 




cd^^d^hF -hdUh 


(A2) 


Then we have 


(A3) 


a = 1 + —^ 2 ^ + 2 T 3 n , 

b = —1 — 2^2^ — 2F3n , 

c = 1 — 2Fin — 2 ^ 2 ^ , 

d = — 1 + 2^1 □ + 2 “^^^ ’ 

/ = 2Fin + F2n-i-2F3n. 


It is easy to see that the form factors a, 6, c, d, f satisfy 
the identities 


a + 6 = 0, c + d = 0, 6 + c + / = 0. 


(A4) 


The equations of motion, obtained from (A2), are 


a(n)n V + K^)dAduh^^ + d^K^) 

+ c{B){r]^^u^p^ah^'" + dpduh) + r]pj,d{n)\Jh (A5) 
+ f{n)n-^d^d,dpd,hp'^ = - 2 KaT^,. 


Here 


2 dSyiatter 

^9, 


(A6) 


flV' 


Let us remember that is a metric in D-dimensional 
Minkowski spacetime and partial derivatives and the □ 
operator are written in Cartesian coordinates in this 
space. Let us emphasize that the number of independent 
arbitrary functions of the □ operator, as well as the form 
of the equations, is the same as in the four-dimensional 
case. However, the dimensional gravitational coupling 
constant depends on the number of dimensions. 

We also show in Sec. H that the form of the equations 


Let write the gravitational potential 'ipd in terms of the 
modes in momentum space 


^dix) = j 


{ 2 ttY 


A'k.x 


V’(k). 


(Bl) 


Here k = A:* is the d-dimensional vector of momentum. 
Similarly we have 

d-^k 


D. 


i{x,x') = J 


( 27 r) 


ilc^x-x') 


Z)(k). 


(B2) 


From (3.15) and ( |3.6[ ) one can derive 
tp{k) = -KdmD{k), 


D{k) = 


1 


= D{k^), 


k = \k\. 


Using the spherical symmetry of the system we get 
d^k = dkdO k^~^ sin^“^ 0 Ad-2^ 


Ad -2 — 2 


7r(c^-i)/2 

u¥y’ 


where Ad -2 is the area of a unit sphere 8 “^ In spherical 
coordinates 


ki{x^ — x'^) = kr cos{d), r=\x — x'\. 
Then the Green function reads 

^(d-l)/2 ^oo ^d-3 


(B3) 


Dd(x,x') = 2 


V)/2 r 

?) Jo 


(27r)‘'a(-fc2) 

X [ dd sin '^“2 0 e*''’'‘=°"(®). 

Jo 

Integration over 0 gives the expression for the Green func¬ 
tion Dd{x, x') in terms of an integral from Bessel function 






- — 1 
o ^ 


ka{—k^) \27rrJ 


(B4) 

Both the potential Yd and the Green function Dd{x,x') 
depend only on the distance r between points. Change of 
the integration variables leads to the following equivalent 
forms {z = kr) 


Dd{r) = 


f 


dz - 


d-4 

Z 2 


( 2 ttYG r^-‘^ Jq a{— z‘^ jr‘^) " 


JdYz). (B5) 


and [t] = z^ jY) 

(B6) 
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